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Abstract 

The partition function of the random energy model at inverse temperature f3 is 
a sum of random exponentials Zi^{l3) = Y^k=i^^'P^P^^-^k)-, where Xi,X2,.-- 
are independent real standard normal random variables (= random energies), and 
n = logA^. We study the large N limit of the partition function viewed as an 
analytic function of the complex variable j3. We identify the asymptotic structure 
of complex zeros of the partition function confirming and extending predictions 
made in the theoretical physics literature. We prove limit theorems for the random 
partition function at complex /?, both on the logarithmic scale and on the level of 
limiting distributions. Our results cover also the case of the sums of independent 
identically distributed random exponentials with any given correlations between 
the real and imaginary parts of the random exponent. 

1. Introduction 

Since the pioneering work of Lee and Yang [3?l[3n], much attention in the statisti- 
cal physics literature has been paid to studying partition functions of various models 
at complex values of parameters such as complex inverse temperature or complex 
external magnetic field, see, e.g., OE]. These studies are sometimes referred to as 
the Lee- Yang program. The motivation here is to identify the mechanisms causing 
phase transitions of the model under study. These transitions manifest themselves 
in the analyticity breaking of the logarithm of the partition function which, in turn, 
is related to the complex zeros of the partition function. Phase transitions are thus 



^Research supported in part by the European Commission (Marie Curie fellowship, project 
PIEF-GA-2009-251200). 



2 



associated with the accumulation points of the complex zeros of the partition func- 
tion on the real axis, in the large system limit. In this respect, complex-valued 
parameters provide a clean framework for identification of phase transitions. 

The main emphasis of Lee- Yang program was on the classical lattice models 
of statistical mechanics. In this work, we focus on the simplest model of a spin 
glass [m [7j: the random energy model (REM) introduced by Derrida [iTJ [12]. 
Let X, Xi, X2, ■ ■ ■ be independent real standard normal random variables. The 
partition function of REM at inverse temperature /3 is defined by 



(1.1) 2N{f3) = Yl 



N 



e 

k=l 



Here, TV is a large integer, and we use the notation n = log N . For real inverse 
temperature /? > the asymptotic behavior of Ziq{j3) as iV — 00 (or equivalently, 
as n — > 00) has been extensively studied in the literature; see [T71 [371 [HI [21 [Z]- 
Specifically, the limiting log-partition function is given by the formula 

(1.2) ,OT:=J.„ ilog|Z„CT| J'+f ' °^''|^- 

A^-i-oo n V 2p, p > v2. 



Convergence (1.2 1 holds both a.s. and in L^, q > 1\ see [T7l[37l[7|. On the level 
of fluctuations, it has been shown in [S] that upon suitable rescaling, the random 
variable Zm{(3) becomes asymptotically Gaussian, for (3 < V2/2, whereas, for /3 > 
it has hmiting stable non-Gaussian distribution, as N ^ 00. 




Figure 1. Complex zeros of Zn in the large A'^ limit. There are three 
phases: Bi (white, no zeros), B2 (light gray, density of zeros is of order 
1), -B3 (dark gray, density of zeros is of order n). On the boundary of Bi 
the linear density of zeros is of order n. The plot shows also the contour 
lines (gray curves and lines) of the log-partition function p. 
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Using heuristic arguments, Derrida |13j studied the REM at complex inverse 
temperature /3 — a + ir. He derived the following logarithmic asymptotics extend- 



ing (1.2) to the complex plane: 

ri + i(a2-r2), /3eSi, 

(1.3) p(/?) lim -log|ZAr(/3)| - <^ V2|a|, /3 e S2, 

where Bi, B2, B3 are three subsets of the complex plane (see Figure [T]) defined by 

(1.4) Bi =. C\B2 U B3, 

(1.5) B2 = {f5eR^:2a^>l,\a\ + \T\>V2}, 

(1.6) B3 = {/? e M^. 2cr2 < l,cr2 > 1}. 

Here, A denotes the closure of the set A. Note that the limiting log-partition 
function p is continuous. 



To derive (1.3), Derrida ^T5] used an approach which can be roughly described 



as follows. Instead of ZAr(/3), one can consider the truncated sum 

N 
k=l 

Indeed, with high probability it holds that Z'^{/3) — Zn{P) since the order of the 
maximum term among . . . , |Xjv| is y/2n and the existence of an outlier satis- 
fying \Xk\ > \^^2n has probability converging to 0. Note, however, that although 
Z^{j3) and ZN[f5) are close in probability, their expectations (and standard devia- 
tions) may be very different from each other, at least for some values of /?. Derrida 
derived an asymptotic formula for the expectation of 2^^(/3), as A'^ — > 00, using the 
saddle-point method. Two cases are possible: the expectation is dominated by the 
energies inside the interval (— \/2n, \/2n) (equivalently, the contribution of the 
saddle point dominates the expectation), or by the energies located near one of the 
boundary points ±-\/2n. He also obtained two similar cases for the standard devia- 
tion of Z'^{/3). Comparing the resulting four formulas, Derrida discovered the three 
phases Bi, B2, B^. The arguments of Derrida [13] are not fully rigorous, although 
it should be emphasized that he did not use the replica method or other standard 
non-rigorous spin glass method. In the present paper, we will make the argument 
of Derrida rigorous and refine his results by deriving distributional limit theorems 
for the fluctuations of Zm{P) (and for the fluctuations in some more general mod- 



els, see Section 2.3 ) at complex /3. An essential feature of the REM at complex 
temperature is the possibility of canceling of terms in Z]\r{/3) due to the presence 
of complex amplitudes. It is for this reason that some standard techniques of rigor- 
ous spin glass theory [43' like the concentration inequalities or the second-moment 
method do not (or do not always) lead to the desired result. These difficulties will 
be discussed in more detail in Section [3.31 



Based on his formula (1.3) for the limiting log-partition function, Derrida [T5j 
computed the asymptotic distribution of zeros of in the complex plane. His 
predictions were in a good agreement with the numerical simulations of Moukarzel 
and Parga [32]. Derrida observed that since Zjv(/3) is an analytic function of /3, 
its empirical distribution of zeros (a measure assigning to every zero weight 1, 
with multiplicities) is given by ^Alog|Zjv|, where A = + -§-2 denotes the 
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Laplace operator in the /3-plane. Taking the large N limit, Derrida obtained the 
formula for the asymptotic distribution of zeros of Zjy. Since the function p 

is harmonic in Bi and B2, Derrida predicted that "there should be no zeros (or at 
least the density of zeros vanishes) in phases Bi and -62"- In phase -B3, "the density 
of zeros is uniform" and is asymptotic to Also, since the normal derivative of p 
has a jump on the boundary of Bi, but has no jump on the boundary between Bi 
and B^ "the boundaries between phases Bi and B2, and between phases Bi and 
B3 are lines of zeros whereas the separation between phases B2 and B3 is not". 
The argument of Derrida involves interchanging the Laplace operator and the large 
N limit. In the present paper we justify Derrida's approach rigorously and derive 
further results on the distribution of zeros of Z^. Namely, we relate the zeros of 
Zpf to the zeros of two random analytic functions: a Gaussian analytic function G 
(in phase B^), and a zeta-function Cp associated to the Poisson process (in phase 
-62)- Also, we will clarify the local structure of the mysterious "lines of zeros" on 
the boundary of Bi. 

For the partition function of REM, considered as a function of a complex external 
magnetic field, a non-rigorous analysis similar to that of Derrida 13^ has been 
carried out by Moukarzel and Parga [331 El] • For directed polymers with complex 
weights on a tree, which is another related model, the logarithmic asymptotics ( 1.3 ) 
has been derived in [T3]; see also [H [22] • Recently, Takahashi [?2] and Obuchi 
and Takahashi ^Hj studied the complex zeros in the generalized REM and other 
spin glass models using the non-rigorous replica method. However, spin glasses at 
complex temperature have not been much studied rigorously in the mathematics 
literature. Our aim is to fill this gap. 

Substantial motivation for the setup with complex random energies comes from 
quantum mechanics. There, the sums of random exponentials with complex-valued 
exponents arise naturally in the models of interference in inhomogeneous media |141 
[15], and in the studies of the quantum Monte Carlo method [16]. 

The sum of random exponentials Z^ is a natural random analytic function ex- 
hibiting, despite of its simple form, a rather non-trivial behavior. We hope that 
the methods developed to study this function can be applied to other random ana- 
lytic functions, for example to random polynomials or random Taylor series. For a 
recent work in this direction, we refer to [Ml [27]. Also, Z^ can be interpreted as a 
(normalized) characteristic function of the i.i.d. normal sample Xi, . . . ,Xjv. This 
connection will be discussed in Section [2^ 

The paper is organized as follows. After introducing some notation in Sect ion |2?l] 



we state our results on zeros and fiuctuations in Sections 2.2 and 2.3 respectively. 
Proofs can be found in Sections [3] and [4] In Section |2.4[ we discuss possible exten- 
sions and open problems related to our results. 



2. Statement of results 

2.1. Notation. We will write the complex inverse temperature /3 in the form (3 ~ 
a + ir, where a, r € M. We use the notation n = logiV, where iV is a large integer 
and the logarithm is natural. Note that in the physics literature on the REM, it is 
customary to take the logarithm at basis 2. Replacing /3 by /3 / -\/log 2 in our results 
we can easily switch to the physics notation. 

We denote by A'r(0, s^) the real Gaussian distribution with mean zero and vari- 
ance > 0. By Nc{0,s'^), we denote the complex Gaussian distribution with 
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density 

w.r.t. the Lebesgue measure on C. Note that Z ^ Nc{0, s^) iS Z = X + iY, where 
X,Y Nu{0,s'^/2) are independent. In this case, EZ ^ and E|Zp = 1. Real 
or complex normal distribution is referred to as standard if s = 1. The standard 
normal distribution function is denoted by $. 

Convergence in probability and weak (distributional) convergence will be denoted 

P w , 

by — > and — >, respectively. Let C be a generic positive constant whose value will 
change at different occurrences. 

2.2. Results on zeros. Let Zn be the partition function of the REM defined as 



in (1.1 1. Note the distributional equalities 



(2.1) Zjv(/3) = (-/?), Zn{I3) = ZN{f3) 



Due to (2.1), it is often enough to consider the case cr, r > 0. The next result 
describes the global structure of complex zeros of Z^, as — ^ oo. Let S3 be the 
Lebesgue measure restricted to B3. Also, let S13 be the one-dimensional length 
measure on the boundary between Bi and B3 (which consists of two circular arcs). 
Finally, let S12 be a measure having the density V2\t\ with respect to the one- 
dimensional length measure restricted to the boundary between Bi and B2 (which 
consists of four line segments). Define a measure S — 2S3 -I- S12 -|- S13. 

Theorem 2.1. For every continuous function / : C — > M with compact support, 
(2.2) - J2 /(/^) .7^ ^ / /(/3)S(d/3). 

Remark 2.2. As a consequence, the random measure assigning a weight 1/n to 
each zero of Z^ converges weakly to the deterministic measure ^'^S. The limit 



measure S is related to the limiting log-partition function p, see (1.3), by the 
formula S = Ap, in accordance with [T3]. Here, A is the Laplace operator which 
should be understood in the distributional sense. The point-wise Laplacian of p is 
easily seen to be 2S3. However, in the distributional Laplacian there are additional 
terms which come from the jumps of the normal derivative of p along the boundaries 
Bi n B2 and Bi CiBs. On the boundary B2 n B3 the jump turns out to be 0. that 
p can be viewed as the two-dimensional electrostatic potential generated by the 
charge distribution S. 



Theorem 2.1 makes the last formula in [13] rigorous. In the next theorems, we 
will investigate more fine properties of the zeros of -Ejy. We start by describing 
the local structure of zeros of Z^ in a neighborhood of area 1/n of a fixed point 
/3o G B3. Let {G{t) : t e C} be a Gaussian random analytic function ^3S] given by 

00 , 

(2.3) G(t) = Va^, 



k=0 ^ "'^ 

where ^Qi^It ■ ■ are independent standard complex Gaussian random variables. The 
complex zeros of G form a remarkable point process which has intensity 1 /tt and 
is translation invariant. Up to rescaling, this is the only translation invariant zero 
set of a Gaussian analytic function; see [HJ Section 2.5]. This and related zero sets 
have been much studied; see the monograph [5T]. 
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Theorem 2.3. Let (5o e B3, be fixed. For every continuous function / : C — )• M 
with compact support, 

^ — ' N^oo ^ — ^ 

^eC: Zn(0)=o ;9eC:G(^)=o 

Remark 2.4. Equivalently, the point process consisting of the points \/n{l3 — /3o), 
where /3 is a zero of Z^, converges weakly to the point process of zeros of G. 

Derrida [13 predicted that the set Bi should be free of zeros. As we will see 
below, it is not true that the number of zeros in Bi converges to in probability 
since with non-vanishing probability there exist zeros very close to the boundary 
of Bi. However, a slightly weaker statement is true. 

Theorem 2.5. Let K be a compact subset of Bi. Then, there exists £ > depending 
on K such that 

V[Zn{P) = 0, for some (3 £ K] ^ 0{N-^), N ^ 00. 

As a consequence, the number of zeros of in K converges to in probability. 
It is natural to conjecture that the convergence holds a.s. The number e, as provided 
by the proof of Theorem |2.5[ converges to as the distance between K and the 
boundary of Bi gets smaller. So, the a.s. convergence does not follow from a Borel- 
Cantelli argument. 

Consider now the zeros of Zn in the set i?2- We will show that in the limit 
as A'^ — ^ 00 the zeros of Z^ in B2 look like the zeros of certain random analytic 
function ^p. This function may be viewed as a zeta- function associated to the 
Poisson process. It is defined as follows. Let Pi < P2 < • ■ • be the arrival times 
of a unit intensity homogeneous Poisson process on the positive half-line. That is, 
Pk = £1 + ■ • ■ + £fe, where ei,S2, ■ ■ ■ are i.i.d. standard exponential random variables, 
i.e., P[ek > t] = e~*, t > 0. For T > 1, define the random process 

(2.4) (p{(3;T) = A^lp.e[o,T] - t'^dt, /3 e C. 

k=i "fc 

Theorem 2.6. With probability 1, the sequence Cp(/3; T) converges as T —> 00 to a 
limiting function denoted by Cp(/3). The convergence is uniform on compact subsets 
of the half-plane {/3 G C: Re/3 > 1/2}. 

Corollary 2.7. With probability 1, the Poisson process zeta-function 

00 

(2-5) Cp(/3) -Y.^ 

fe=l '^k 

defined originally for Re/3 > 1, admits a meromorphic continuation to the domain 
Re/3 > 1/2. The function C,p{P) — Cp{P) ~ ^giii ^-S- analytic in this domain. 

The next theorem describes the limiting structure of zeros of Z^ in B2. The 
form of the process C,p appearing there is not surprising and can be explained as 
follows. In phase B2 the process Zm is dominated by the extremal order statistics 
of the sample Xi, . . . , Xjv- These form a Poisson point process in the large N limit, 
see, e.g., [521 Corollary 4.19(i)], and Qp is some functional of this process. 
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Theorem 2.8. Let f : B2 ^ R be a continuous function with compact support. Let 

C,p^ and Cp'' be two independent copies of Qp. Then, 

E E m+ E /(-/^)- 

C<.^'(/3/y2)=0 C?'(/3/y2)=0 



Theorem 2.8 tells us that the zeros of Z^y in the domain a > 1/v 2, \a\ + \t\ > v2 
(which constitutes one half of B2) have approximately the same law as the zeros of 
Cp, as N OD. Let us stress that the approximation breaks down in the triangle 
a > I/V2, |<t| + |t| < v^- Although the function Cp is well-defined and may have 
zeros there, the function Zff has, with high probability, no zeros in any compact 
subset of the triangle by Theorem |2.5| 

Next we state some properties of the function (p. Let (3 > 1/2 be real. For 
/? ^ 1, the random variable Cp(/3) is stable with index and skewness parameter 



1. In fact, (2.4) is just the series representation of this random variable; see [HI 
Theorem 1.4.5]. For /3 = 1, the random variable Cp(1) (which is the residue of (p 
at 1) is 1-stable with skewness 1. For general complex /3, we have the following 
stability property. 

Proposition 2.9. If Cp \ ■ ■ ■ i Cp^'' independent copies of C.p, then we have the 
following distributional equality of stochastic processes: 

C^p'^ + ... + Ci'^ ^k^Cp- 

To see this, observe that the union of k independent unit intensity Poisson pro- 
cesses has the same law as a single unit intensity Poisson process scaled by the factor 
1/fc. As a corollary, the distribution of the random vector (Re Cp(/3), Im Cp(/3)) be- 
longs to the family of operator stable laws; see [5T] . 



Proposition 2.10. Fix r G M. As a l 1/2, we have 
V2ct - 1Cp(ct + ir) - 



As a corollary, there is a.s. no meromorphic continuation of C,p beyond the line 
(7 = 1/2. Using the same method of proof it can be shown that for every different 
Ti,T2 > the random variables \/2a — 1 (p{a+iTj), j = 1,2, become asymptotically 
independent as cr J, 1/2. Thus, the function ^p looks like a nai've white noise near 
the line cr = 1/2. The intensity of complex zeros of (p at f3 can be computed by 
the formula g{f3) = ^AElog |Cp(/3)|, where A is the Laplace operator; see dU 



Section 2.4]. Proposition 2.10 suggests that g(/3) ~ ^ (2cr~i)^ as cr J, 1/2. In 
particular, every point of the line cr = 1/2 should be an accumulation point for the 
zeros of ^p with probability 1. 

Let us look locally at the zeros of Zjv near some /^o = cro + ^^o on one of the 
boundaries Bi n B^ or Bi O B2. We will show that in both cases the zeros form 
approximately an arithmetic sequence. The structure of the measures S13 and S12 



in Theorem 2.1 suggests that the distances between the consequent zeros should 
behave like — in the first case and like in the second case. The next theorems 

n \To\n 

show that this is indeed true. First, we analyze the boundary Bi D B3. 
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Theorem 2.11. Let /3o = ctq + itq be such that CTq + Tq = 1 and ctq < 1/2. There 
exist a complex-valued random variable ^ and a bounded real sequence (5jv such that 
for every continuous function / : C — > M ivith compact support, 

/3gC:2n(/3)=0 ^ ^ ^ fceZ 

Remark 2.12. In other words, the zeros of Zjy near /3q are given by the formula 

^ ^ / 27ri/s + ^ + i(5Ar \ /1\ 
/3 = /3ofl + ^ 

As we will see in the proof, the random variable Re^ takes negative values with 
positive probability. It follows that the probability that has a zero in Bi does 
not go to as — oo. 



The boundary Bi n B2 consists of 4 line segments . By symmetry (2.1), it suffices 
to consider one of them. 

Theorem 2.13. Let Pq — ao + iTQ be such that ao > 1/V2, tq > andao+Tg = \f2. 
There exist a complex-valued random variable r/ and a complex sequence d^ — 
0(log n) such that for every continuous function / : C — > M with compact support, 

/3gC:2jv(/3)=0 fcSZ ^ ^^'0 

Remark 2.14. In other words, the zeros of Zn near /3q are given by the formula 

^ ^ 2x, 1 / 2nik \ / 1 \ , „ 

^ = /3o +6" 3 - + rf^rj +o l^-j , keZ. 

2.3. Results on fluctuations. We state our results on fluctuations for a general- 
ization of ( |1.1[ ) which we call complex random energy model. This model involves 
complex phases and allows for a dependence between the energies and the phases. 
Let {X, Y), {Xi, Yi), . . . be i.i.d. zero- mean bivariate Gaussian random vectors with 

Var Xk = Var n = 1, Corr(Xfc, Yk) = p. 



Here, — 1 < p < 1 is fixed. Recall ( 1.1 1 and consider the following partition function: 

N 



(2.6) Zjv(/3) = ^e^("^^+*"^'=\ (3 = (a,T) e 



k=l 



For T = 0, this is the REM of Derrida [12 at real inverse temperature a. For p — 1, 
we obtain the REM at the complex inverse temperature j3 — a -\- ir considered 
above; see ( |1.1[ ). For p = 0, the model is a REM with independent complex 
phases considered in [14] . Note also that the substitutions (/3,p) {—(3,p) and 
{f3,p) (/3, —p) leave the distribution of Zn{(3) unchanged. 
Recall (1.2). Define the log-partition function as 

(2.7) pjv(/3) = -log|Zjv(/?)|, f3={a,T)e 



1)2 



n 



Theorem 2.15. For every /3 e M^, the limit 
(2.8) p(/3) lim p^iP) 

N^oo 
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exists in probability and in L'^ , q>l, and is explicitly given as 

■l + i(a2-r2), (3eB,, 
(2.9) p{f3) ^ { V2\a\, (3eB2, 



Note that the hmit in ( |2.9[ ) does not depend on p. However, we will see below 
that the fluctuations of Zn{/3) do depend on p. The next theorem shows that 
Zn{/3) satisfies the central limit theorem in the domain < 1/2. 

Theorem 2.16. If < 1/2 and t ^ 0, then 

(2.10) ZnW)~N ^ 

Remark 2.17. If < 1/2 and r = 0, then the limiting distribution is real normal, 
as was shown in 8J. 

Remark 2.18. If in addition to cr^ < 1/2 we have cr^ > 1, then N'^+hi'^^-''^) = 
o{N2+'^ ) and, hence, the theorem simplifies to 

(2.11) 1^ A^c(0,l). 



Eq. (2.11) explains the difference between phases Bi and B^: in phase Bi the 
expectation of Zn{/3) is of larger order than the mean square deviation, in phase 
B^ otherwise. 

In the boundary case = 1/2, the limiting distribution is normal, but it has 
truncated variance. 



Theorem 2.19. If a"^ = 1/2 and t ^ 0, then 



iVc(0,l/2). 



N 

Next, we describe the fluctuations of Zisi{f3) in the domain > 1/2. Due to 



(2.1 ), it is not a restriction of generality to assume that cr > 0. Let ^at be a sequence 
such that V^bNG^"^'^ ~ as iV — > cx). We can take 

(2.12) b^ = V2^- 



2V2n 

Theorem 2.20. Let a > 1/^2, r ^ 0, and \p\ <l. Then, 
Z^(/3)-jVE[e^Wry)lx<b.] ^ ^ 

where Sa denotes a complex isotropic a-stable random variable with a characteristic 
function of the form E[e*^'=(^°^)] = q-^^o'^MA" ^ z(eC 

Remark 2.21. If cr > l/\/2 and r = 0, then the limiting distribution is real totally 
skewed a-stable; see |8j. If cr > l/\/2 and p = 1 (resp., p = —1), then it follows 
from Theorem 14.81 below that 



ZM{P)-NneP^^lx<bA J- ? ( P\ ( ? ( 
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Remark 2.22. We will compute asymptotically the trmicated expectation on the 



2:201 



left-hand side of (2.13|) in Section 3.2 below. We will obtain that under the as- 
sumptions of Theorem 



(2.15) ^^^5 iia+\T\>V2, 

(2.16) 2N{fi)-N . iia+\r\<V2. 

Similarly, if cr > l/\/2, but p = 1, then we have 

For p — —1, we have to replace /? by /3. 

2.4. Discussion, extensions and open questions. The results on fluctuations 

are closely related, at least on the heuristic level, to the results on the zeros of 



Zn- In Section 2.3 we claimed that regardless of the value of /3 7^ we can find 
normalizing constants mAr(/3) G C, vpf{P) > such that 

ZAr(/3) - TOjv(/3) ^ 

— ^ '-'^^ 

for some non-degenerate random variable Z{(3). It turns out that in phase Bi the 
sequence mN{P) is of larger order than vn{P)^ which suggests that there should 
be no zeros in this phase. In phases B2 and i^a, the sequence mis[il3) is of smaller 
order than fAr(/3), which does not rule out the possibility of zeros in these phases. 
One way to guess the density of zeros in phases B2 and B3 is to look more closely at 
the correlations of the process Zn. In phase ^3, it can be seen from Theorem |4.6| 
below that Zn{13i) and Zn{I32) become asymptotically decorrelated if the distance 
between /3i and (32 is of order larger than l/y/n. This suggests that the distances 
between the close zeros in phase i?3 should be of order 1/-\Aj and hence, the density 
of zeros should be of order n. Similarly, in phase B2 the variables Zn{(3i) and 



Zn{P2) remain non-trivially correlated at distances of order 1 by Theorem 4.8 
below, which suggests that the density of zeros in this phase should be of order 1. 

An additional motivation for studying Zn comes from its connection to the em- 
pirical characteristic function. Given an i.i.d. standard normal sample Xi, . . . , X^r, 
the empirical characteristic function is defined by Cn{(3) = X^feLi e^^'^'' . We have 
Zn{P) = CN{—i\/n(3). The limit behavior of the stochastic process {cjv(/3) : /3 € K} 
without rescaling (3 by the factor ^/n has been much studied; see, e.g., [TFl llOj. 
There has been also interest in the behavior of Rn = inf{/3 > 0: Re cat (/3) = 0}, 
the first real zero of Recjv; see [SUl [22]. In particular, it has been shown in [201 
Corollary 4.5] that, for aU t e M, 

lim P[R% -n<2t] = $(-\/2e-*). 

N—^oo 

Hence, the first real zero of ReZjv(/3) restricted to /3 G iM. is located near i with 
high probability. This is exactly the point where the imaginary axis meets the set 
B3. 
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It is possible to extend or strengthen our results in several directions. The 
statements of Theorem |2.1| and Theorem |2.15| should hold almost surely, although 
it seems difficult to prove this. Several authors considered models involving sums 
of random exponentials generalizing the REM; see [HI US IHj • They analyze the 
case of real /3 only. We believe that our results (both on zeros and on fluctuations) 
can be extended, with appropriate modifications, to these models. 

3. Proofs of the results on fluctuations 

3.1. Truncated exponential moments. We will often need estimates for the 
truncated exponential moments of the normal distribution. In the next lemmas, we 
denote by X a real standard normal random variable. Let $(2;) = —j^ f e~~da: 
be the distribution function oi X. It is well-known that 

1 + 0(1) .2 

(3.1) $(z) = ^=Me-T, 

The normal distribution function $ can be extended as an analytic function to the 
entire complex plane. We need an extension of (3.1) to the complex case. 



Lemma 3.1. Fix some e > 0. The following holds as \z\ — > 00, z e C; 
(3.2) $(z) 



-^e-V, z/|argz|>f +£, 



In particular, $(z) — > 1 z/ |z| — > 00 and \ argz| < j — e. 

Remark 3.2. We take the principal value of the argument, ranging in (~7r,7r] and 

having a jump on the negative half-axis. In the domain ^ -I- e < | argz| < ^ — e 

both asymptotics in (3.2) can be apphed. To see that they give the same result, 

,2 

note that I -e~T" I 00 there. 



Proof of Lemm a\3l\ For the first case of ([3^, see P, Eq. 7.1.23 on p. 298]. The 
second case of ( |3.2[ ) follows from the identity $(2) = 1 — $(—2:). □ 

In the next lemmas, we record several simple facts on the truncated exponential 

2 

moments which we will often use later. Note that E[e"'"^] = for all w € C. 

Lemma 3.3. Lei w e C, a e M. Then, E[e"'^lx<Q] = e'#'I>(a - w). 
Proof. For w G M, we have 

E[e"'-^ljf<a] = / e^'^-Tdz = ^e^ / e'^^dz = e^$(a - w). 
V27r J -00 V27r J-00 

For w G C, this holds by analytic continuation. □ 
Lemma 3.4. Let a G M. The following estimates hold. 

(1) Ifw>a, thenE[e'"^lx<a]<e°-'"~'^ ■ 

(2) Ifw<a, thenE[e'"^lx>a] < e'""-^. 

Proof. Consider the case w > a. By Lemma 3.3 E[e™'^lx<a] = e^^{a — w). 



Using the inequality $(2) < e '2 valid for 2 < 0, we obtain the statement of 
case (1). Case (2) can be reduced to case (1) by the substitution {X,w,a) H> 
{-X,-w,-a). □ 
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Lemma 3.5. Let F{n) — Ele^^^ 1 x ^^a{n)] ' where w = u + iv£C, 72>0, 
and a{n) is a real-valued function with lim„_5.oo a{n) = a. The following hold, as 
n ^ oo: 

(3.3) 



F{n) i \/27rn(to-a) 



^n{a{n)w-^a («))^ if U + \v\ > O, 

if u + \v\ < a. 



If w eR and a{n) — w + + o(^), for some c G M, then 
(3.4) i^(n) - $(c)e5"''", n -> oo. 



Remark 3.6. The second hne in ( |3.3[ ) can be generahzed to the foUowing formula 
valid in the case u — \v\ < a: 



(3.5) 



F{n) 



1 + "(1) n(a(n)u)-ia2(n)) 



v27rn(w — a) 

Proof of Lemma \3.5\ Let 2;(n) = •y/na(n) — Wy/n. By Lemma 3.3 

'^"$(\/na(n) - w\/"^) 



we have 



(3.6) 



F{n) 



<i>(z(n)). 



Note that 2:(ri) ^ {a — u — iv)^/n, as n — ?■ oo. 

Case 1. If M + |ti| > a, then | arg z{n)\ > j + e, for some e > 0, and all sufBciently 



large n. Applying the first line of (3.2 1, we arrive at the first line of (3.3) 



Case 2. If u — \v\ < a, then | argz(n)| < ^ — s, for some £ > 0, and all sufficiently 



large n. Applying the second hne of we get (3.5). 



Case 2a. If even the stronger condition ■u+|i'| < a holds, then ^{u'^—v'^) > au 



1^2 



and hence, the first term in (3.5) asymptotically dominates the second one. We 



obtain the second line of (3.3). 

Case 3. If a = w € M and a(n) ~ w ^ — ^ + o(^), for some c G M, then 



lim„_i.oo z{n) — c and we arrive at (3.4) 



□ 

Lemma 3.7. // {X,Y) is a real Gaussian vector with standard marginals and 
correlation p, then, for s, a G M, 



X<a\ 



X<a\ 



In particular, E[e^('"-^+*^'*')] 



Proof. We have a distributional equality {X,Y) = {X,pX + — p-^W), where 
{X, W) are independent standard normal real random variables. It follows that 



Efe 



s{<jX + iTY)-: 



'■X<a\ 



X<a\ 



_ ^Us(a + iTp)X+iSTy^l-p^W J 



where we have used that E[e*^] = e*'/^ and that W and X are independent. □ 



3.2. Proof of Theorems 2.16 , 2.19, 2.20, The main tool to prove the results 



on the fluctuations is the summation theory of triangular arrays of random vec- 
tors; see [19] and [31]. The following theorem can be found in [19l §25] in the 
one-dimensional setting and in [30] or in [31] Theorem 3.2.2] in the d-dimensional 
setting. Denote by | - 1 the Euclidean norm and by (•, •) the Euchdean scalar product. 
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Theorem 3.8. For every N (z N, let a series Wi_n, ■ • • , Wn.n of independent ran- 
dom d-dimensional vectors he given. Assume that, for some locally finite measure v 
on M'^\{0}, and some positive semidefinite matrix S, the following conditions hold: 

(1) liniAT^oo Ef=i ^Wk,N £ B]= v{B), for every Borel set B C M''\{0} such 
that v{dB) = 0. 

(2) The following limits exist: 

N N 



E = limlimsup^ Var[VFfe,Arl|H/^ „|<e] 



limliminf V Var[iyfc,Arl|vi/, „|<e] 

k=l 



Then, the random vector Sn '■= SfeLi(^fc,iv ~ IE[Wfcjvl|vi^^. „|<ii:]) converges, as 
N ^ (X), to an infinitely divisible random vector S whose characteristic function is 
given by the Levy-Khintchine representation 

I- i{t,s)l\s\<il)v{ds), te 



logE[e*<*'^>] 



-(t,^t) 



Here, R > is any number such that v does not charge the set {s e 



= R}. 



Proof of Theorem 2.16 For k ~ 1, . . 



, N , define 



/n{aXk+iTYk) 



Let Wn be a random variable liaving tli e sa me law as the Wk,N^&- Note that 
iVE[VFAr] = 7v(i-'^'-^'+2''"^'')/2 by Lemma 
show that 



3.7 



N 



■ N) 



k=l 



JV->oo 



To prove the theorem, we need to 
iVc(0,l). 



The proof is based on the two-dimensional Lindeberg central limit theorem. We 
consider Wk^N as an M^-valued random vector (Re VFfc_Ar, Im VF^ jv). Let Ejy be the 
covariance matrix of this vector. First, we show that 



lim NT,M = 



A/2 



N^oo " V 1/2. 



(3.7) 
We have 

(3.8) NE[{ReWN)^ + (ImVKjv)^] = A^E[|V7ivn = N-'^''^E[e^''^^] = 1. 
Also, we have iVE[W^^] = jy-^r^+'Harp Lemma 3.7 Since we assume that t ^ 0, 
this implies that limAr_>.oo -/VE[W^] — 0. By taking real and imaginary parts, we 
obtain that 

(3.9) lim NE[{ReWNf - {ImWNf 



N 



lim NE[{ReWN)ilmWN)] = 0. 



Combining (3.8) and (3.91, we get 
(3.10) 



lim NE[{ReWNy]= lim iVE[(ImW0v)2] = 1/2. 



Also, by Lemma 3.7 we have 
(3.11) lim 



N- 



NE\Wn] = lim N-^^" +^--2i^rp)/2 ^ 



It follows from pTLOf , ( |3TT| ) that (|37| holds. Fix an arbitrary e > 0. We 

complete the proof of the theorem by verifying the Lindeberg condition 



(3.12) 



N 



lim NE[\Wn-EW, 



]=0. 
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Assume first that cr 7^ 0, say a > 0. Write = cr- 



2a 



(7+2^ > 2(7 by the assumption cr^ < 1/2. Hence, by Part^of Lemma 



Then, lini 



3.4 



(3.13) 



N 



lim ArE[|W^Ar|2l 



\Wn\>b_ 



~ lim 



2(T n|grg2(TyriJ\r J 



we have 
0. 



This also trivially holds for cr = 0. Together with ( |3.11[ ), ( |3.13| ) imphes ( |3.12[ ). □ 

-Proo/ 0/ Theorem [2.19[ Without loss of generality, let a = l/\/2. For fc = 1, . . . , iV, 
define 



N 



Let Wn be a random variable with the same distribution as the Wu^n^S- To prove 
the theorem, we need to verify that 



N 

E 

k=l 



{Wu,N-^Wk,N) iVc(0,l/2). 

N^oo 



As we will see in equation (|3.14|) below, the Lindeberg condition (|3.12|) is not 
are going tc 
By Lemma 



3.5 



satisfied. We are going t o ap ply Theorem |3.8| instead. Fix e > and let = 



(3.14) 



equation (3.4) with c = 0, we have 



1 



X<i^/na j\ 



1/2. 



If IpI ^ 1, then by Lemma |377| and (3.14) 



(3.15) lim iVE[W^^l|vi.„|<e] 



< lim e-2"(i-''')^'iV-iE[e 



2nX^ 



N 



0. 



The result of (3.15) continues to hold for \p\ = 1 since in this case, Lemma 3.7 and 
Lemma 3.5 (first part of (|3.3[)) yield, as — > 00, 



NE[W'^ll 



By Remark |3.6| we have, as A^ — cx) 



X<^/na 



o{N- 



-1 ra(y2ajv-3a^)^ 



0. 



E[e- 



fiiaX -1 



. V2 



1 2 ^ 



It follows that 
(3.16) 



lim iVE[|W^^|l|H.„|>e] - Jim^E[e^-^l^>^,„] = 0. 



We consider Wj^ as an M^-valued random vector (Re M^jv, Im VFat). It follows 
from (3.14), (3.15), (3.16) that the covariance matrix Ejv '■— Var[Vl^jvl|Wjv|<£] sat- 
isfies 



(3.17) 



lim A^Ejv = 



A/4 



It follows from (3.16) that limjv-i-oo A^P[|WAr| > e] = 0. Therefore, the conditions of 
Theorem 3.8 are satisfied with v = Q and E given by the right-hand side of (3.17). 
Applying Theorem |3.8[ we obtain the required statement. □ 



Prooj of Theorem \2J^ RecaU that a = G (0, 2). For fc = 1, . . . , A^, define 

random variables 
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Let Wn be a random variable having the same law as the Wk.N's. We will verify 
the conditions of Theorem |3.8| To verify the first condition, fix < ri < r2, 
< di < 02 < in and consider the set 



B = {z e C: ri < \z\ < r2,0i < arg z < 6*2}. 



We will show that 
(3.18) 

Define a set 
We have 



lim NV\Wn e B] 



2n 



A 



N 



u R 



C 



P[Wn eB]= p[e'^v^(-Y-''«) (n, r2), r e An] 

= /" ' P [y e Ajv I fT\/^(X - ^at) = logr] fN{r)dr. 

Jri 

Here, fN{r) is the density of the log-normal random variable e\^'^(^^^"): 



(3.19) /jv(r) = 



1 



■ exp ■ 



1 / log r 



+ h 



N 



N 



N ->■ 00, 



where the asymptotic equivalence holds uniformly in r G [?'i,7'2]. To prove ( |3.19 1, 
recaU that V27r6jve''"/2 N and bj^ ^ V2n. Conditionally on a^/n{X — 6jv) = 
logr, the random variable Y is normal with mean ^j-n = + ^n) and variance 

v/1 The variance is strictly positive by the assumption \p\ 7^ 1. It follows 

easily that 

lim ¥[Y e An \ a^{X - 5w) = logr] = 



Bringing everything together, we arrive at (3.18). So, the first condition of Theo- 
rem [3]8] holds with 

a dxdy 



2tt r2+" ' 

To verify the second condition of Theorem |3. 8| with S = 0, it suffices to show that 



(3.20) 



lim limsnvNE[\WN\^l\w«\<e 



]=0. 



Condition |W^Ar| < e is equivalent to X < qn, where un = bN + ^ V2n. By 



Lemma 



3.5 (first case of (3.3)) with A = ^/n, w = 2a, we have 



where we have again used that -\/27r5Are''"/^ ^ N. We obtain that 



hm 7VE[|W^jv|'l|i^„|<e] 



lim iVe~2av/Sb„]£jg2^v^^j^^^^] ^fj^2-V2/a 



Recalling that 2 > v2/cr, we arrive at (3.20). By Theorem 

N 



3.8 



]) .7^ 
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where the Hmiting random vector Sa is infinitely divisible with a characteristic 
function given by 

, dxdy 



V'(z) :=logE[e*<^°'^>] 



a 
2^ 



(e 



i{u,z 



1 — i(u, z) 1 



\u\<l 



\2+a ■■ 



z e 



Here, u = x + iy and {u, z) — Re(uz). Clearly, i^{z) depends on \z\ only and satisfies 
ip{Xz) = A"'0(z) for every A > 0. It follows that ■ip{z) — const • □ 



Proof of Remark\KM_We prove pls] ) and ( pie] ) first. Let a > 1/^/2, r 7^ 0, and 
\p\ < 1. By Lemma 3.7 we have 



(3.21) 



rriN 



7VE[e^^(''^+^^^)lx<6„] = A^l-5^'(l-p')E[e^^('^+'^'')^l^ 



<bj\ 



Write w = (J + irp. Recall from (2.12) that 
(3.22) V2^6Are''"/2 ~ TV, 



N 



2n, iV 



Applying Lemma 3.5 we obtain 

(3.23) E[e^(^+*^'')^lx<6„] ^ <^ ^^^i^^-' 

In the case a + |rp| = v^, we applied Remark 3.6 
in (3.5 1 dominates the second one. 



and noted that the first term 



Proof of ( [2I5I ). Assume that cr+|r I > ^2. If even the stronger condition (t+|tp| > 
\/2 is satisfied, then it follows from ( 3.21[ ) and the first line of (3.23) that 

(3.24) Irrijvl - Ce-^v^fciv-ir^Ci-p')" ^ oie"-^^"). 

The last step follows from t 7^ and \p\ < 1. If a + |t | > V2 but cr + \tp\ < V2, 
then it follows from (3.21) and the second line of (3.23) that 

(3.25) |m^| -iVi+5('^'-"') ^o(( 



The last step follows from the inequality l + ^(c r^— r^) < ^/2a. It follows from (3.24) 
and (3.251 that we can rewrite Theorem 2.20 in the form ( |2.15 1. □ 



Proof of \2.1Q\ . Assume that cr + |t| < ^2. Then, a - \tp\ < ^2 and it follows 
from (3.21) and Remark 3.6 that 

mjv = Ari-5^'(i-f')(e3'"'" + 0(A^-ie'"^^")) = N^+h(-^''^'}+i''^p + o(e''^''"). 



The last step follows from t 7^ and \p\ < 1. Hence, we can rewrite Theorem 2.20 
in the form (12. 161). □ 



We now proceed to the proof of (2.17) and (2.18). Let a > l/\/2 and p = I- 
Our starting point is (2.14). 



Proof of (2.17). Assume that ct + |t| > and a 7^ \/2. By Lemma 
of (3.3 1, we have 

V2 



3.5 



first line 



rriN 



N 

/2^(/3- \/2)' 



2 '^N 



where we have used (3.22). Recall that 
(3.26) Cf 



72 
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It follows that we can rewrite (2.14) in the form (2.17) 



□ 



3.6 



Proof of ( |2.18[ ). Assume that a + \t\ < V2. By Remark 

:= 7VE[e^v/^^l^<,„] - N'+'^^' ^ + o(l) 



rriN 



It follows that we can rewrite (2.14) in the form (2.18) 



□ 
□ 



3.3. Proof of Theorem |2.15[ We will deduce the stochastic convergence of the 
log-partition function p7v(/3) = ^ log |Zjv(/3)| from the weak convergence of Zn{/3). 
This will be done via Lemma [3^ stated below. One may ask whether there exists a 
more direct way to prove the convergence of pm{P)- A standard method to handle 
such questions for real /? is to use the Gaussian concentration inequality; see |431 
Theorem 1.3.4]. To apply it we need to verify that the function pAr(/3) is Lipschitz in 
the variables Xi, . . . , Xj^. This is easy to do in the real setting, but if /3 is complex, 
we have Zn{13) = for some non-empty set of tuples Xi,. . . ,Xpf. Thus, pn{P) 
is not even finite, so that the Lipschitz property does not hold. The possibility of 
having infinite pjv(/3) is not just a technical difficulty, especially in view of the fact 
that the zeros of Zn become dense in B-^ in the large N limit. 



Lemma 3.9. Let Zi, Z2, ■ ■ . he random variables with values in C and let tojv G 
C, ujv G C\{0} be sequences of normalizing constants such that 

(3.27) Z^-^N ^ ^ 

Vpf N^oo 

The following two statements hold: 

log|Zw| P 



(1) // |uAr| = oQuinI) and |to7v| 00 as iV 00, then " 1. 

(2) // |mjv| — 0{\vi\[\), \vn\ —>■ c» as iV —> 00 and Z has no atoms, then 

logl^ivl I 
log 1 1' JV I jv^^, 

Proof of (1). Fix e > 0. For sufficiently large N , we have |tojv| > 1 and, hence, 



i-e<;^<i+. 

log|TOAr| 



> 





Zn - mN 


1 


\mN\ 






<2 






VN 




\vn\ . 



The right-hand side converges to 1 by our assumptions. 
Proof of (2). Fix e > 0. For sufficiently large N, 



□ 



log|t;jv| 



> 1 



\vn\ 



> \vn\' 



< p 




Zn - mN 














VN 





The right hand-side converges to by our assumptions. Consider now 



log I WAT I 



< 1 



\vn\ 



< \vnI 



niN 



VN 



■mN 

VN 



< IvnI' 



Assume that there is (5 > such that the right-hand side is greater than 5 for 
infinitely many A^'s. Recall that tun/vn is bounded. Taking a subsequence, we 
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may assume that —misf/vN converges to some a G C. Recall that \vn\ — > oo. But 
then, for every ry > 0, 

V[\Z - a| < 77] > hmsupP 

JV->oo 

This contradicts the assumption that Z has no atoms. □ 



Zn - mN 



VN 



< 



> S. 



Proof of Theorem 2.15[ Convergence in probability. Let p(/3) be defined by (2.9 1. 
Note that p is a continuous function. We are going to prove that for every /3 g C, 
limTv-foo PAf(/^) — p{l^) in probability. We may assume that r 7^ since otherwise 



the result is known from equation ( 1.2 1; see [T71I32- follows from Theorems 2.16 
2.19[ |2.20[ and Remark [2^ that condition ( |3.27[ ) is satisfied with Zn = 2n{I3) 
and an appropriate choice of normalizing sequences tunt^n- We will now verify 
that Lemma |3.9| is applicable. 

Case 1. Let < 1/2 . In this case, tun and vn are given by Theorems 



2.18 



Namely, \mN\ = N^+hi'^^-^^) and vn = N^' 



2.16 



2.19 



see also Remark 

Case la. If in addition to cr^ < 1/2 we have cr^+r^ < 1, then l + i(cr^ — r^) > 
and we obtain \vn\ = o{\mN\)- 

Case lb. If in addition to cr^ < 1/2 we have cr^ + r^ > 1, then l + |((T^-r2) < |+cr^ 
and we obtain \mN\ = 0{\vn\)- 

Case 2. Let cr^ > 1/2 and, withou t restriction of generality, a > 1/^/2. Then, m^r 

gcr^/nbjv g^j^jj ^jjg formula for niN 



and Vn are given by Remark 2.22 Namely, \vn\ 
depends on whether cr + |r| > \/2 or not. 

Case 2a. If cr > 1/V2 and ct+ |t| > \/2, then rriN = 0, see ( |2.15| ) and ( |2.17[ ). Thus, 
o(|wAr|) is satisfied. 



TUN 



Case 2b. If cr > l/^/2 and cr + |r| < V2, then |77ijv| = N^+H'^"'^^^); see (pie]) 



-hi-' 



> 



\/2cr it follows that |wAr| = o(|tojv|) 



and ( 2.18 1. From the inequality 1- 
Case 2c. If cr > 1/^2 and cr + |t| %/2, then 1 + ^(cr^ - t^) = ^/2cr. However, 
y/2n —00 by (2.12 1, we still have \vn\ — o(\mN\). 



since 



fnb 



N 



To summarize, the normalizing constants itln and vn satisfy the first condition 
of Lemma 3.9 if G -Bi or /3 belongs to one of four (open) line segments on the 



bounda ry o f Bi and i?2. Otherwise, tojv and vn satisfy the second condition of 



Lemma 



3.9 



Note that we need also to verify that the ran dom v ariable (p{f3/\/2) 
has no atoms if cr > l/y/2. This will be done in Lemma 3.10 below. Applying 
Lemma 3.9 we obtain that pn{P) ^ Pif^) in probability. □ 



Lemma 3.10. If u > 1/2, then the random variable C,p{P) ^0,^ ""^o atoms in C. 

Proof. For a random variable Y with values in C, let Q{Y) — supj^gj,P[y — y] 
be the weight of the maximal atom of Y . Note that Q is a special case of the 
concentration function; see [551 p. 22]. For independent random variables Yi and 
Y2^ the convolution formula implies that 



(3.28) 



Q{Yi+Y2) < max(Q(ri),Q(r2)). 



Also, Q{Y + c) = Q{Y) for every c e C. Let Pi < P2 < ■ • • be the points of a unit 
intensity Poisson point process on the positive half- line. For T > 1, we write 



(3.29) 



CpW)^CpW;T) + Ri/3;T), 



19 



where R{f3;T) is a rest term and C,p{f3]T) is defined as in (2.4), that is 



(3.30) 



p 

k=l 



Note that (p{(3;T) and R{(3;T) are independent random variables since R{f3;T) 
depends only on those points of the Poisson process which are in the interval (T, oo). 

We will show that Q{Cp{P; T)) < e"^ for every T > 1. By ( |3.28[ ) and ( |3.29[ ) this 
implies that Q{Cp{f3)) = 0, which is the desired result. Consider the random events 
Am{T) — {J2T=i lpfcG[o,T] — G Nq. That is, Am{T) occurs if there are 

exactly m points of the Poisson point process in [0,T]. Note that P[Ao(T)] = e~'^ 
and (p{f3;T) is constant on the event Ao{T). Let z E C By the total probability 
formula, 



(3.31) 



■[Cp{(3;T) ^ z] < 



J2nCp{P;T)^z\A„ 

m—l 



, P„i have the same dis- 



Conditionally on where m e N, the points Pi 
tribution as the increasingly reordered independent random variables Ui 
distributed uniformly on [0,r]. It follows that, for every m € N, 



(3.32) 



P[Cp(/3;T) = z\A„,] <Qij2Uk^]< QiUi^) = 0, 



^fc=i 



where the last line follows from the fact that the random variable Ui ^ has no 
atoms. It follows from ( [33T| ) and ([332]) that P[Cp(^;T) = z] < e"^ for every 
2; e C, T > 1. This implies that Q{Cp{/3; T)) < e^^ and completes the proof. □ 

Proof of Theorem 2.15\ Convergence in i'. We are going to show that pn{P) 
p{(3) in L'^, where g > 1 is fixed. From the fact that pn{I3) — >■ p{f}) in probability, 
and since p(/3) > 0, for every /3 G C, we conclude that, for every C > p(/3). 



lim Pw(/3)lo<p„(/3)<c+i =P{I3) in L''. 
For every w € M, we have 

V[pn{(3) >u]< e-™E|2:Ar(/3)| < e-""iVE[e'"^^] = e"(<^-"), 
where C = 1 + (t^/2. From this, we conclude that 

C30 

E [|m(/3)rip„(/3)>c+i] = 5^E[|p^(/?)|«lc+fe<p„(/3)<c+fc+i] 

00 

fe=i 

which converges to 0, as — > 00. To complete the proof, we need to show that 
(3.33) lim E[|p^(/3)rip„(^)<o] =0. 

A* — >-oo 

The problem is to bound the probability of small values of Zn{(3), where the log- 
arithm has a singularity and |p7v(/5)| becomes large. This is non-trivial because of 
the presence of complex amplitudes in the definition of Zff{(3); see (2.6). We have 
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to show that there is not much canceUation among the terms in (2.6). Fix a small 
£ > 0. Clearly, 

(3.34) E [|pAr(/3)ri_e.^<p„(^,)<o] < {e<J^r. 



To prove (3.33), we would like to estimate from above the probability P[|Zjv(/3)| < r] 
for < r < e-^" Recall that Zjv(/3) is a sum of N independent copies of the 
random variable Q\/ri.{'^x+iTY) ^ Unfortunately, the distribution of the latter random 
variable does not possess nice regularity properties. For example, in the most 
interesting case p = 1 it has no density. This is why we need a smoothing argument. 
Denote by Br{t) the disc of radius r centered at t G C. Fix a large A> 1. We will 
show that uniformly in i G C, l/A < < A, n > (2A)2, and < r < e"^''''", 

(3.35) P[e^('"-^+*^^) e B^(t)] < Crw. 

This inequality is stated in a form which will be needed later in the proof of Theo- 
rem HHJ 

Let \t\ > and r > 1/(2^). The argument argt of a complex number t is 
considered to have values in the circle T = K/27rZ. Let P : M — !> T be the canonical 
projection. Denote by Ir{t) the sector {z G C: |arg2; — argt| < 2y^}, where we 
take the geodesic distance between the arguments. A simple geometric argument 
shows that the disc Br{t) is contained in the sector Ir{t). The density of the random 
variable P{Ty/nY) at e T is given by 



By considering the right-hand side as a Riemann sum and recalling that r > 1/{2A), 
we see that the density converges to l/(27r) uniformly in G T as iV oo. We 
have 

which implies (3.35). 

Let now \t\ < ^/r. Then, recalling that logr < —ea'^n, we obtain 

log r 



consider tl 
variable e"^^ is given by 



X < 



_(logr)^ ^ 
< C 16<T^ri < r 16 . 



It remains to consider the case t > y'r, \<t\ > 1/(2A). The density of the random 



1 _(logx)^ 

5(2;) = ^== — c 2.^^ , x>0. 
V Zirnax 

It attains its maximum at xq = The maximum is equal to ^(a^o) — 

i_ga^n/2_ ^ < (27rn)cr2e-'"'". Then, 

Let r > (27rn)f72e"'^^", which, together with \a\ > 1/(2^) and n > (2^)^, implies 
that r > Using the unimodality of the density g and the inequality t — r>r, 

we get 

p[gv^(^x+»ry) ^ < p[t - r < e"^^ <t + r]< 2rg{r) < c"^^^ < . 
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The last inequality follows from r < e 



This completes the proof of (3.351 



Now we are in position to complete the proof of (3.33). Let Ur be a random 



variable distributed uniformly on the disc i?.r(0) and independent of all variables 



considered previously. It follows from (3.35 1 that the density of the random variable 
^^{aX+iTY) _|_ jj^ bounded from above by Cr-^+f^/^^). Hence, the density of 
Zn{P) + Ur is bounded by the same term C'r^^+'^Z^"^. With the notation r = e"*"'" 
it follows that, for every k > ecr^. 



< 



\ZN{I3) + Ur\ < 2r] < Cr'to. 



P[pA'(/3)<-fc]-P[|ZAr(/3)|<e- 
From this, we obtain that 

n\PN(fi)np.{P)e[-k-i.-k]] < C{k + l)'?e-*. 
Taking the sum over all k — ea^ + /, Z = 0, 1, . . ., we get 



E[K(/3)rip„(0)< 



< Ce- 
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l = \ 
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Recalling (3.34), we arrive at (3.33) 



□ 



Remark 3.11. As a byproduct of the proof, we have the following statement. For 
every A > 0, there is a constant C = C{A) such that E|pjv(/3)| < C, for all 
l/A <\P\ < A and sufficiently large N. 



4. Proofs of the results on zeros 

4.1. Convergence of random analytic functions. In this section, we collect 
some lemmas on weak convergence of stochastic processes whose sample paths are 
analytic functions. As we will see, the analyticity assumption simplifies the things 
considerably. For a metric space Af , denote by C{M) the space of complex- valued 
continuous functions on AI endowed with the topology of uniform convergence on 
compact sets. Let 13 C C be a simply connected domain. 

Lemma 4.1. Let {[/(<): t G 13} he a random analytic Junction defined on D. Let 
T d D be a closed differ entiable contour and let I'C be a compact subset located 
strictly inside F. Then, for every p € Nq, there is a constant C ~ Cp{K,T) such 
that 



E 



sup|[/(p)(t)| 
teK 



<C i ¥.\U{w)\\dw\ 



r 



Proof By the Cauchy formula, \U^P\t)\ <C §^ \U{w)\\dwl for all t G K. Take the 
suprcmum over t Cz K and then the expectation. □ 

It is easy to check that a sequence of stochastic processes with paths in C (D) is 
tight (resp., weakly convergent) if and only if it is tight (resp., weakly convergent) 
in C{K), for every compact set K C D. 

Lemma 4.2. Let Ui, U2, ■ ■ ■ be random analytic functions on D. Assume that there 
is a continuous function /: Z3 —> M such that E|t/jv(i)| < fit), for all t € D, and 
all N Cz N . Then, the sequence Ujq is tight on C{D). 
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Proof. Let K C D he a compact set. Let F be a contour enclosing K and located 
inside D. By Lemma ITT] 



E 



sup I [/at (i) I 



<C (j) fiw)\dw\, E 



sup I C/^ (01 



<C f f{w)\dw\ 



By standard arguments, this implies that the sequence Un is tight on C{K). □ 

Lemma 4.3. Let U, Ui, U2, . . . be random analytic functions on D such that Un 
converges as N ^ 00 to U weakly on C{D) and P[C/ = 0] = 0. Then, for every 
continuous function f : D ^ R with compact support, we have 

zGC: Un{z)=0 zeC: (7(z)=0 

Remark 4.4. Equivalently, the zero set oi Un, considered as a point process on D, 
converges weakly to the zero set of U. 

Proof. Let H be the closed linear subspace of C{D) consisting of all analytic func- 
tions. Consider a functional 5* : iJ — > M mapping an analytic function ip which 
is not identically to /(•^)' where the sum is over all zeros of tp. Define also 
^"(0) — 0. It is an easy consequence of Rouche's theorem that 5* is continuous on 
iJ\{0}. Note that iJ\{0} is a set of full measure with respect to the law of U. 
Recall that Un ^ U weakly on H. By the continuous mapping theorem |39i §3.5], 
'^'{Un) converges in distribution to 'i'{U). This proves the lemma. □ 



4.2. Proof of Theorem |2.1[ A standard approximation argument shows that we 
can assume that / is infinitely differentiable. Let A be the Lebesgue measure on 
C. In his computation of the limiting density of zeros, Derrida [T3] used the fact 
that ^Alog |Zjv| (where A is the Laplacian interpreted in the distributional sense) 
gives the measure counting the zeros of Zn- That is, 

(4.1) E = ^ /log|Z^(/3)|A/(/3)A(d/3). 

/3eC: 2n(^)=o 

A proof of ( |4.1[) c an be found in 21, Section 2.4]. Recall that p{f3) has been defined 
in Theorem 12. 151 We have 



(4.2) / p(/3)A/(/3)A(d/3) = / /(/3)S(d/3). 

Jc Jc 

Indeed, Green's identity gives 

p(/?)A/(/3)A(d/3) = 1^ Ap{i3)fmm+£^ {p((^^^^ - f(f^^^) 

Here, n is the unit normal to the boundary of Bi pointing outwards Bi and is 
the corresponding directional derivative. The first term on the right-hand side is 
equal to 2 f {/3)'E.^{d/3) for i = 3 and to for i = 1, 2. Adding Green's identities 
for i = 1,2, 3, noting that ^ has no jumps and computing the jumps of ^ on the 



B, 



boundaries between the different B^'s, we arrive at (4.2) 



Recah that pn{I3) = Mog |ZAr(/3)|. From (|4.1| and ((42]) we conclude that 
Theorem |2.1| is equivalent to 

pjv(/3)A/(/3)A(d/3) A / p(/3)A/(/3)A(d/3). 

N^oo 
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We will show that this holds even in . By Fubini's theorem, it suffices to show 
that 



(4.3) 



lini 

AT-foo 



E|m(/3)-p(/3)||A/(/3)|A(d/3)=0. 



We know from Theorem 2.15 that limAr^oo E|PAf (/?) — p(/3)| = 0, for every /3 G C. 
To complete the proof, we need to interchange the limit and the integral. We may 
represent / as a sum of two functions, the first one vanishing on |/3| < 1/4 and 
the second one vanishing outside |/3| < 1/2. Since the contribution of the second 
function to (2.2) vanishes by Theorem 2.5 we may assume that / vanishes on 
1/3 1 < 1/4. With this assumption, the use of the dominated convergence theorem is 
justified by Remark |3.11[ □ 



4.3. Proof of Theorem 2.5 , The idea of the proof is to show that the fluctuations 
of Zn{I3) around its expectation are of smaller order than the expectation, in phase 
Bi. We don't rely on the expression for the limiting log-partition function p. Let 
r be a differentiable contour enclosing the set K and located inside Bi. We have 



P[Zjv(/3) = 0, for some P E K] < 



sup 

fieK 



ZNW)-EZNm 



< E sup 
<C <j)E 



Zw(/3)-EZw(/3) 



> 1 



EZn{P) 



\dw\ 



where the last step is by Lemma 



4.1 



Note that \EZn{(3)\ - N^+ki'^'-^') . To 
complete the proof, we need to show that there exist £ > and C > depending 
on F such that, for every /3 eV, N eN, 

(4.4) E|Zjv(/3) -EZjv(/3)| < CN^-'+'^^^^-^^l 

Since F C we can choose e > so small that F C B[{e) U B"{e), where 
B[{e) = {/3 e C: CT^ + < 1 - 2e}, 
B'lie) = {/3 e C : (|a| - V2)^ - > 2e, 1/2 < < 2}. 

We have 

E\ZnW) ~EZnW)\^ = A^E|e^^-^ -Ee'^^^p < iVEe^'"^^ = N^^^"'. 
If /3 e B[{e), then it follows that 

E \Zn{/3) - EZn{P)\ < A^H'^' < ^l-e+i(^'-r^)^ 



This impHes (4.4). Assume now that /3 G B'{{e) and a > 0. For k — 1, . 
define random variables 



Uk,N — C -'^Xfc<V2Tr' 

By Part [T] of Lemma |3.4[ we have 



'Xfc>V2n' 



N 



(4.5) E 



k=l 



k.N 



EUk. 



< NElUi 



N 



Ne 



-2^2(7 



X<V2n 



]<1. 
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Similarly, by Part [2] of Lemma |3.4| 



(4.6) E 



N 



k=l 



Combining (4.5| and (4.6), we obtain E\Zn{P) ~EZn{P)\ < 3e' 



Since P e 



B'{{e), this implies the required estimate (4.4). □ 
As a by-product, we obtain a proof of the formula for the limiting log-partition 
function in phase Bi which is simpler than the proof given in Section |3.3[ We will 
use only the information about the first two truncated moments of Zjsf(/3). 



Proposition 4.5. For (3 £ Bi, we have limAr_j.oo 
probability. 



PnW) = 1 + ^(^2 - r2) in 



Proof. We have shown in (4.4) that, for every /3 e there exist C > and e > 
depending on /3 such that for all N E N, 



(4.7) 



E 



- 1 



< CN- 



Let p(/3) = 1 + i(cr2 _ r2). Note that ^ log |EZAr(/3)| = p(/3) for all N e n. It 
follows that for every i5 > and all sufficiently large iV, 



(4.8) P[bw(/3)-p(/3)| > 5] 



log 



Zn{P) 



EZn{P) 



> nS 



< 



ZnW) 



EZnW) 



1 



> 



In the last step, we have used that |log|2;|| > nS implies that |z — 1| > 1/2, for 
large n. By (4.7) and the Markov inequality we can estimate the right-hand side 
of (4.8) by 2CN~'^ , which implies that pn{P) converges to p{(3) in probability. □ 

4.4. Proof of Theorems |2.3| and 2.1H Recall that G is the Gaussian analytic 
function defined in (2.3). Theorem 2.3 will be deduced from the following result. 

Theorem 4.6. Fix some [3o = ao + iTQ with ctq < 1/2 and tq ^ 0. Define a random 
process {GNit) : i e C} by 



(4.9) 



GN{t) := 



-■N 



N 



Then, the process Gn converges weakly, as N oo, to the process e * /^G(t) on 
C(C). 



Proof. For k = 1, 



, N, define a random process {Wk.N{t) : t G C} by 

Wk Nit) = iV-V2g(/3o\A^+t)Xfc-(<Tov^+t)^^ 



Then, Gjv(t) = J2k=i(^k,N{t) — ^Wk,Nit))- First, we show that the convergence 
stated in Theorem 14.61 holds in the sense of finite-dimensional distributions. Take 
ti, . . . ,td e C. Write = (Wk.Niti), ■ . ■ , Wk,N{td))- We need to prove that 



N 



(4.10) 



fc=l 



k.N 



EWk,N) ^ (e- 

Af-i-oo 



-tl/2 



<G{ti 



,e''"-/^G{ta)). 
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Let Wn be a process having the same law as the W^fe.Ar's and define Wjv = 
{WN(ti), . . . , W]y{tci)). A straightforward computation shows that for aU s e C, 



(4.11) 
(4.12) 

Also, we have 
(4.13) 



NE[WN{t)WN{s)] 



lim N\E[WN{t)WN{s)] \ = 0. 



lim 

AT— foo 



N\E[WN{m 



lim e"5('^o+^o)"+0(v^) 



0. 



Note that by (2.3 1, 



E[e-* /2G(t)e--V2G(s)] 



We see that the covariance matrix of the left-hand side of (4.10) converges to the 
covariance matrix of the right-hand side of (4.101 if we view both sides as 2d- 
dimensional real random vectors. To complete the proof of (4.10), we need to 
verify the Lindeberg condition: for every e > 0, 



(4.14) 



N 



hm 7VE[|W^|2l|w„|>e] =0. 



For ^ = 1, . . . , d, let A/ be the random event |M^Ar(ti)| > |W^jv(ij)| for all j 
On Ai, we have [W nI"^ < d\WN{ti)\'^ . It follows that 



iVE[|W 



N\ 



MWjv|>e 



] < d"^NE 

1=1 



\WN{tl)\H^W.^U)\> 



^0, 



where the last step is by the same argument as in (3.13). This completes the proof 
of the finite-dimensional convergence stated in (4.10). The tightness follows from 
Lemma 4.2 which can be applied since 

E\GN{t)\ < vnc^W] < vmiwNim = ep-*)'. 

The last equality follows from (4.11). □ 

Proof of Theorem \2.3\ If (5q e B^, then the expectation term in the definition of 
Gjv, see (4.9), can be ignored: we have limjv_j.oo |G'Ar(t) — C/jv(t)| = uniformly on 
compact sets, where 



UN{t) = N- 



l3o 



t 



It follows from Theorem 



4.6 



that Um converges to the process e *^/^G(t) weakly 

□ 



on C(C). Applying Lemma 4.3 we obtain the statement of Theorem 2.3 



Proof of Theorem \2.11\ Let 5n be a bounded sequence such that na^TQ — SM G 27rZ. 



Taking t = /3o^y|^ in Theorem 



4.6 we obtain that weakly on C(C), 

s + iSn ' 



N- 



G(0). 



Doing elementary transformations, we arrive at 
2^(/3o(l+^)) 



G(0). 



The zeros of the right-hand side are located at s = 2nik + S^, k E Z, where ^ = 

□ 



log(— G(0)). The proof is completed by applying Lemma 4.3 
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4.5. Proof of Theorems [2T6t [2T8| and \2A3[ 



Proof of Theorem \2.6\ Fix a compact set K contained in the half-plane a > 1/2. 
Define random C(i4r)-valued elements Sk{l3) = si(/?) + . . . + Sfe(/3), where 

Note that si, S2, . . . are independent. By the properties of the Poisson process, 



(4.15) 



E[sfc(/3)] = 0, ^E[|sfc(^)p] = / r^-dt < 

k=l 



Thus, as long as ct > 1/2, the sequence {5'fe(/3)}fcgN is an L^-bounded martingale. 
Hence, Sk{f3) converges a.s. to a limiting random variable denoted by S{f3). We need 



to show that the convergence is uniform a.s. It follows from (4.15) and Lemma 4.2 
that the sequence Sk, k G N, is tight on C{K). Hence, Sk converges weakly on C{K) 
to the process 5*. By the Ito-Nisio theorem [24^, this implies that 5*^ converges to 
S a.s. as a random element of C{K). This proves the theorem. □ 

Proof of Theorem \2.8\ Let us first describe the idea. Consider the case a > l/-\/2. 
Arrange the values Xi, . . . , in an increasing order, obtaining the order statistics 
Xi-N < • ■ • < Xn;n- It turns out that the main contribution to the sum Zn{(3) = 
J2k=i e^^-^'' comes from the upper order statistics Xjs-k-.N, where fc = 0, 1, . . .. 
Their joint limiting distribution is well-known in the extreme-value theory, see [391 
Corollary 4.19(i)], and will be recalled now. Denote by M the space of locally finite 
counting measures on E = K U {-|-cx)}. We endow M with the (Polish) topology of 
vague convergence. A point process on M is a random element with values in M. 
Let Pi , P2 , • ■ ■ be the arrivals of the unit intensity Poisson process on the positive 



half-line. Define the sequence bi^ as in (2.12), that is 



, /TT- log(47rn) 



2V2n 



Proposition 4.7. The point process ttn '■— ^iV^i^k ~ ^n)) converges as 

N ^ 00 to the point process tToo = Sfe°=i ^(~(log^/£)/V2) weakly on M. 

Utilizing this result, we will show that it is possible to approximate Z]y{f3) (after 
appropriate normalization) by C,p{P/y/2) in the half-plane a > 1/^/2- Consider 
now the case a < — 1/-\/2. This time, the main contribution to the sum Zn{(3) 
comes from the lower order statistics Xk-.N, k — 1,2, .. . Their joint limiting dis- 
tribution is the same as for the upper order statistics, only the sign should be re- 
versed. Moreover, it is known that the upper and the lower order statistics become 
asymptotically independent as iV — 00; see or [321 Cor. 5.28]. Thus, in the 
half-plane a < —1/^/2 it is possible to approximate Z^iP) by an independent copy 
of (^p{— (5 / \/2) . In the rest of the proof we make this idea rigorous. For simplicity 
of notation we restrict ourselves to the half-plane Z? = {/3eC: cr> l/\/2}. 

Theorem 4.8. The following convergence holds weakly on C{D): 
^^^^^ r^rL 1,71 
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The proof consists of two lemmas. Take A > and write ^nIP) — Cjv(/5) 
e^*(/3)+A^(/3), where 

N 
k=l 

e^(/3) = iVE \^e^V^^^>'^f>^)i^^_ 

N 



= E 



E 



^/3yn(Xfc-6]v)-i 



k=l 



Lemma 4.9. Lei CpIS') ^6 defined as in (2.4|. Then, the following convergence 
holds weakly on C{D): 



^^{P)-eUP) Cp 

N^oo 



Proof. Recall that by Proposition 4.7 the point process ttjv converges to the point 
process tToo weakly on M. Consider a functional M — C{D) which maps a 
locally finite comiting measure p — J^ieiHUi) G M to the function 'i>{p){f3) = 
^^gj e^^' lj,.>_^, where f3 £ D. Here, / is at most countable index set. If p 
charges the point +oo, define ^'(p), say, as 0. The functional 4' is continuous on 
the set of all p e M not charging the points —A and +oo, which is a set of full 
measure with respect to the law of tToo. It follows from the continuous mapping 
theorem [33 §3.5] that = '^/{ttn) converges weakly on C{D) to '^{tToo)- Note 
that 



fc=l 



We prove the convergence of e^(/3). Using the change of variables ^/n{x — 6jv) = y, 



we obtain 



N 

J2i 



Recalling that -\/27r6Are''«/^ ^ TV and ^ \/2n as iV oo, we obtain that 



lim7v-+oo e^(/3) = t '^/^di, as required. 



□ 



Lemma 4.10. For every compact set K d D, there is C > such that, for all 
sufficiently large N , 



E 



sup \Aim 



< Ce^^"^'^)^/^ 



Proof. Let F be a contour enclosing K and located inside D. First, E[A;^(/3)] = 
by definition. Second, uniformly in /3 e F it holds that 



X<b„ 



\ On 

< Ce(i-^'^)^, 



A 
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where the second step follows from Lemma 3.3 and the last step follows from (3.1 1. 
By Lemma [4J| we have 



E 



sup \A^m 
peK 



<C j E|A^(/3)||d/?| < c'e(i-^'^)^/2. 



The proof is complete. 



□ 



Proof of Theorem 4-8 By Theorem 2.6 we have the weak convergence 



Cp 



V2 



;e^^ 



A— s-oo 



V2 



Together with Lemmas |4.9| and |4. 10[ this implies Theorem |4.8| by a standard argu- 
ment; see for example [261 Lemma 6.7]. □ 



The proof of Theorem 
yields 



2.8 



can be completed as follows. For a > Lemma 



3.5 



f V2 

lim A^e-^^''"E[e''^^lx<fc„] = { 



oo, 



if \a\ + |r| > A 
if M + \t\ < v/2. 



The first equality holds uniformly on compact subsets of B2 ■ By Theorem |4.8[ the 
process e~'^^''"Zjv(/3) converges to Qp(j3/y/2) weakly on the space of continuous 
functions on the set B2 H {cr > 1/^/2}. Similarly, on the space of continuous func- 
tions on n {cr < — 1/\/2} the same process converges w eakly to an independent 
copy of Cp{— 13/^/2). By Lemma 4.3 this implies Theorem 



2.8 



□ 



Proof of Theorem 2.13[ Let c?^ be a complex sequence such that 

d'^ -I- /?o^°^^^^ - iTgW e 27riZ and d'^ = O(logn). 
2v 2 

Write /^AT — Bo+ where s G C is a new variable. Note that limAr_j.oo 0n = 

(Pfl — v2)ri 

/3q. Let X ^ iV]R(0, 1). Applying Remark 3.6 and noting that the second term on 



the right-hand side in (3.51 converges to 



"V2 



, we obtain 



lim A^e-''"^''"E[e''"^^ljf<f,„l lim iVe-''"^^"e^ + ^ ^ 



By Theorem 4.8 and (3.26), the following holds weakly on C(C): 



'■TV 



72 



The zeros of the right-hand side are located at s = 27rifc + rj^ k € where 77 
log(— Cp(/3o/'s/2))- Define — d'j^ / {^/2tq) . The theorem follows from Lemma 
after elementary transformations. 



4.3 



□ 
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4.6. Proof of Proposition 2.10, Let r 7^ be fixed. Let S{f3) be a random 
variable defined as in the proof of Theorem 2.6 Take a,b E R. For a > 1/2, 
consider a random variable 

(00 /.fe 1 

^/(P,;a)li<p^<fc-y f{t;a)dt 



where f{t; <t) = ^/a^ + b'^t °' cos(t log t 
and sin 6 = , , , , . 



and €E M is such that cos 9 



We need to show that — 1 i^(o-) converges, as ct J, 1/2, to a centered real 
Gaussian distribution with variance (a^ + h'^)/2. By the properties of the Poisson 
process, the log-characteristic function of Y{<7) is given by 



logEe'^^^" 



/OO 
f{t;a)dt. 



We will compute the second term and show that the first term is negligible. By 
elementary integration, we have 

a2 + 52 1 + cos(2r log t - 20) 



(4.16) 



fit;a)dt. 



a 



2 

2^52 



-dt 



Using the inequalities |e*- 



(4.17) 



zfit-a) 



2 

1 — ix 
1-i 



1 



-2Bi 



Re 



2(7-1 (l-2CT)+2iT^ 

I < \x\^ and |/(<;cr)| < Ct-" , we obtain 



zf{t-a) + -f{t-a)]dt 



< 



C 



3(7-1' 



Bringing (4.16) and (4.17) together and recalling that r 7^ 0, we arrive at 
(4.18) lim logEc''^''2^^^('^) = - ^ (a^ + h^)z^ . 

This proves the result for r 7^ 0. For t = 0, the limit is (4.18) is —a?z'^/2. 
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